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ABSTRACT 

We report on high resolution Echellc spectroscopy of 20 giant stars in the 
Galactic old open clusters NGC 6791 obtained with Hydra at the WIYN tele¬ 
scope. High precision radial velocity allow us to isolate 15 bona fide cluster 
members. From 10 of them we derive a global [M/H]=+0.39±0.05. We there¬ 
fore confirm that NGC 6791 is extremely metal rich, exhibits a few marginally 
sub-solar abundance ratios, and within the resolution of our spectra does not 
show evidences of spread in metal abundance. With these new data we re-derive 
the cluster fundamental parameters suggesting that it is about 8 Gyr old and 
4.3 kpc far from the Sun. The combination of its chemical properties, age, po¬ 
sition, and Galactic orbit hardly makes NGC 6791 a genuine Population 1 open 
cluster. We discuss possible interpretations of the cluster peculiarities suggesting 
that the cluster might be what remains of a much larger system, whose initial 
potential well could have been sufficient to produce high metallicity stars, and 
which has been depopulated by the tidal held of the Galaxy. Alternatively, its 
current properties may be explained by the perturbation of the Galactic bar on 
an object originated well inside the solar ring, where the metal enrichment had 
been very fast. 

Subject headings: open clusters: general — open clusters: (NGC 6791) 


1. Introduction 

NGC 6791 is an extremely interesting and intriguing open cluster. The combination of 
old age, small distance and high metal abundance makes this cluster very attractive, and 
indeed in the last 40 years it has been the target of intensive and numerous studies (Carney 
et al. 2005, and references therein). A large number of optical photometric studies (Stetson 
et al. 2003 and references therein) has been recently complemented by the deep ACS/HST 
investigation by King et al. (2005), and the near IR study by Carney et al. (2005). 

Since the pioneering study of Kinman (1965) it was clear that NGC 6791 is a very old and 
very metal rich cluster. 
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Its age was measured several times by using different sets of isochrones (Carraro et al. 
1999 and references therein), and it is probably confined in the range from 8 to 12 Gyr, 
depending on the cluster precise metal abundance. Taylor (2001, and references therein) 
critically reviewed all the available metallicity estimates, concluding that the [Fe/H] for 
NGC 6791 should probably lie in the range +0.16 to +0.44 dex. This combination of age 
and metallicity is unique in the Milky Way open cluster population, and has been recently 
questioned by Bedin et al. (2005), whose HST study of the White Dwarf cooling sequence 
supports a much younger age. As these authors comment, this age discrepancy may arise 
from defects in the White Dwarf current models, or from the poorly known cluster metal 
abundance. 


Noteworthy, the cluster is also known to harbour a number of sdB/sdO stars (Landsman 
et al. 1998, Buson et al. 2005), which may be explained by a scenario of a high metallicity 
driven wind in the Red Giant Branch phase of the progenitors of these stars or, more simply, 
by the binarity hypothesis (Green et al. 2005) . 

The UV upturn (namely the abrupt rise in the UV continuum emission shortward of 
A ~ 2000A) similar to that typical of any ellipical galaxy (Landsman et al. 1998) and the 
highly eccentric orbit, unusual for a Population I object, contribute to make this cluster even 
more intriguing . 

In an attempt to substantially improve our knowledge of NGC 6791, we carried out a 
spectroscopic campaign to provide radial velocities and accurate metallicities of a statistically 
significant number of stars in the cluster. In fact, current abundance determinations either 
lack sufficient resolution or are restricted to a very small number of stars. 

This new set of abundance estimates coupled with the high quality of existing photometry 
(Stetson et al 2003) allow us to significantly improve on the fundametal parameters of this 
cluster, and better clarify its intriguing nature. 


2. Observations 

The observations were carried out on the night of July 28, 2005 with the Hydra spectro¬ 
graph at the WIYN telescope at Kitt Peak observatory under photometric conditions and 
typical seeing of 1". 1 arcsec. The MOS consists of the Hydra positioner, which in 20 minutes 
can place 89 fibers within the 1-degree diameter focal plane of the telescope to ~ 0.2 arcsec 
precision. This project employed the 3-arcsec diameter red-optimized fiber bundle. The 
fibers feed a bench-mounted spectrograph in a thermally isolated room. With the echcllc 
grating and Bench Spectrograph Camera the system produces a resolution of 20,000 at 6000 
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A . The wavelength coverage of 200 A around the central wavelength of 6000 A provides 
a rich array of narrow absorption lines. We observed 20 RGB/clump stars with 45 min 
exposures, for a grand total of 4.5 hours of actual photon collection time on the same single 
star. The 20 stars where selected from the Stetson et al. (2003) photometric catalog to be 
giant stars and to have the right magnitudes to be observed with the WIYN 3.6m telecope. 
We restricted the sample to giant stars brighter than V ~ 15. The stars are listed in Table 1, 
where first column reports Stetson et al. (2003) numbering and column 2 Kinman (1965) 
numbering (K65). Then coordinates, magnitudes and colors are taken from Stetson et al. 
(2003). The radial velocities and spectral classification have been derived in this paper, 
following Villanova et al. (2004). 


3. Data Reduction 

Images were reduced using IRAF 2 , including bias subtraction, flat-field correction, frame 
combination, extraction of spectral orders, wavelength calibration, sky subtraction and spec¬ 
tral rectification. The single orders were merged into a single spectrum. As an example, we 
show in Figure 2 a portion of the reduced, normalized spectrum for star #11814 where some 
spectral lines are identified. Some spectra have very low S/N, although all the observed stars 
have practically the same magnitude. This could happen for two reasons: the first one is an 
imperfect pointing of the fiber and the second one a possible bad fiber-trasmission. Because 
of this, we could not use five stars for abundance measurements (see below). 


4. Radial Velocities 

Radial velocities (RV) for RGB and clump stars in NGC 6791 have been determined 
many times in the past. Kinman (1965) obtained radial velocities for 19 stars, and spectral 
type for 21. Later, RVs have been measured by Geislcr (1988, 12 stars), Fricl et al. (1989, 9 
stars), Garnavich et al. (1994, 18 stars), Scott et al. (1995, 32 stars) and Friel et al. (2002, 
41 stars), with different resolution and precision. 

We derived here RVs for 20 stars (see Table 1). The radial velocities of the target stars were 
measured using the IRAF FXCOR task, which cross-correlates the object spectrum with a 


2 IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the As¬ 
sociation of Universities for Research in Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 
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template. As template, we used a synthetic spectrum calculated by SPECTRUM (see 5.2 for 
a description of the program) with roughly the same atmospheric parameters and metallicity 
of the observed stars. The final error in the radial velocities was typically less than 0.2 km 
s _1 , and in many cases less than 0.1 km s -1 . These errors are significantly lower than in 
any other previous investigation. This allowed us to clean out held interlopers and isolate 15 
bona fide members. RVs are plotted in Figure 3. Five stars have radial velocites completely 
different from the others, and so were considered non-members, although it possible that 
some of them are binary stars. 

A few of our targets are in common with previous investigations, and we can have an external 
check on our RV measurements. For 11 stars we provide the first estimate of the RV. 

In general, we find that Garnavich et al. (1994) RV estimates (for stars 3003, 3010, 3036, 
2001, 3018 and 2008) are systematically larger than ours, by about 5-7 km s -1 , although, 
given their typical large error (5 to 15 km/sec), these differences cannot be considered statisti¬ 
cally significant. Also Kinman (1965) RVs for the two stars in common with our investigation 
(stars 2001 and 2008) are larger than our estimate. The largest deviation is with respect to 
the RVs by Fricl et al. (2002), and their previous measurements (Scott et al. 1995, Fricl et 
al. 1989). In this case, there are differences in some case exceeding 20 km s^ 1 (stars 3003 
and 2008). Finally, within the errors, we find a good agreement for the two stars (stars 3003 
and 3010) we have in common with Geisler (1988). 

From the RVs of the 15 cluster members in our sample, we obtain a mean radial velocity 
V r = —47.1 ± 0.8 km s -1 , in good agreement, within the errors, with the values obtained by 
the other authors, with the exception of Kinman (1965). The RV dispersion results to be 
ay = 2.2 ± 0.4 km s -1 . 


5. ABUNDANCE ANALYSIS 
5.1. Atmospheric parameters 

The atmospheric parameters were obtained from the photometric BVIc data of Stetson 
et al. (2003). According to Stetson et ah, the most likely reddening and absolute distance 
modulus are E(B-V ) = 0.09 [(E(V-I) = 0.11)] and (m — M) 0 = 12.79. Effective tempera¬ 
tures ( T e ff ) were obtained from the color-T e // relations of Alonso et al. (1999), Sekiguchi 
& Fukugita (2000), and Ramirez & Melendez (2005). The temperatures, obtained from the 
B-V and V-I colors using the quoted relations, are in agreement within 50-100 °K. The 
gravity log(g) was derived from the canonical formula : 

*<) = 4xM|t)-Mi) + /<) 

In this equation, the mass M/M 0 was derived from Straizys & Kurilicne (1981). The hum- 
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nosity L/L 0 was derived from the absolute magnitude My , adopting the distance modulus 
of Stetson et al. (2003). The bolometric correction (BC) was derived from the BC -T e ff 
relation from Alonso et al. (1999). The typical error in log(g) is 0.1 dex. 

Finally, the adopted microturbulence velocity is the mean of the values given by the relation 

(Gratton et al. 1996): 

v t = (1.19 x 10“ 3 )T e// - 0.90 log(g) - 2; 

and the relation (Hondashelt et al. 2000): 

v t = 2.22 — 0.322 log(g). 

The typical error in v t is 0.1 kms^ 1 . 


5.2. Abundance determination 

The resolution (R= 17,000 at 6580A) of our spectra, the high metallicity of the cluster, 
and the low temperature of the target stars cause a lot of blending, and therefore it was 
not possible to measure the equivalent width of the single spectral lines. For this reason, 
the abundances were determinated by comparing the observed spectra with synthetic ones. 
The synthetic spectra were calculated by running SPECTRUM, the Local Thermodynamical 
Equilibrium (LTE) spectral synthesis program freely distributed by Richard O. Gray (see 
Piotto et al. 2005 for the details on our synthetic spectra calculation). Model atmospheres 
were interpolated from the grid of Kurucz (1992) models by using the values of T e // and log(g) 
determined as explained in Section 5.1. We analized stars with T e // > 3900°/l because, 
for lower temperatures, molecular bands were present, creating difficulties for continuum 
determination. We could analyze only ten stars, after rejecting the non-members and the 
stars which resulted to be too faint and too cool. 

First of all, we compared the entire spectrum (range 6400-6760 A) with the synthethic 
one in order to obtain an estimate od the global metallicity [M/H], Then, we analized single 
lines in order to measure abundances of Fe, Ca, Ti, Ba, Al, Ni, and Si. 

A preliminary lines-list was obtained considering all the strongest lines present in our 
spectra, identified using the line-list distributed with SPECTRUM. The final lines-list was 
created from the preliminary one by comparing the observed solar spectrum with a sinthetic 
one, calculated with SPECTRUM for the Sun parameters (T e // = 5777 K, log(g) = 4.44, v t = 
0.8 m/s). The lines in the synthetic spectrum which did not properly match the observed 
ones were rejected. We also checked whether the preliminary line identification was correct 
using the MOORE line database (Moore et al. 1966). Table 3 shows the final list of lines 
we used in our analysis. The resulting metallicities for each star arc listed in Table 4. 
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Due to the radial velocity shit, a few lines (6462.567, 6475.630, 6493.781, 6532.890 A) 
overlapped with telluric lines. We did not consider these lines in the abundances determina¬ 
tion. 


An example of the comparison between synthetic and observed spectra is shown in 
Fig. 4, for the case of star #11814 Finally, using the stellar parameters [colors, T e ff, and 
log(g)] and the absolute calibration of the MK system (Straizys & Kuriliene 1981), for each 
star, we derived the stellar spectral classification (see Villanova et al. 2004 for details), which 
is listed in Col. 9 of Table 1. 


The weighted mean of the [Fe/H] content of the ten members of NGC 6791 analysed in the 
present paper is [Fe/H] = +0.39±0.01 (internal error). Previous investigations reported a 
variety of estimates for the metallicity of NGC 6791. By using low resolution spectroscopy, 
Fricl & Janes (1993) obtained [Fe/H] = +0.19T0.19 from 9 stars, and Friel et al. (2002) 
[Fe/H] = +0.11±0.10 from moderate resolution spectra of 39 stars. Because of the large 
errors, the first estimate is compatible with our one, within one sigma, while the second one 
is off by almost 3a. 


Interestingly enough, our results are in very good agreement with the study by Peterson 
& Green (1998), who derived [Fe/H] = +0.40T0.10 for star #2017, a cool blue horizontal 
branch star, using a resolution very similar to the one used in the present study. 


In conclusion, our results confirm that NGC 6791 is actually a very metal rich cluster. 
Noteworthy, within the errors of our measurements, the metallicities listed in Table 4 do not 
show any significant abundance spread. 


5.3. Abundance ratios 

Abundance ratios constitute a powerful tool to assign a cluster to a stellar population 
(Friel et al. 2003, Carraro et al. 2004, Villanova et al. 2005). 

In Table 5 we list the abundance ratios for the observed stars in NGC 6791. These values do 
not show any particular anomaly. All the abundance ratios are solar scaled, with the only 
exception of [Al/Fe] and [Ba/Fe], which seem to be slighlty under-abundant. Our ratios are 
in good agreement with those provided by Peterson & Green (1998). 



6. Distance and age of NGC 6791 


Our accurate determination of the metal content of NGC 6791 allows a new, more 
reliable estimate of the cluster distance and age. 

To this purpose, in this section, we are going to fit the observed Color Magnitude Diagram 
(CMD) from Stetson et al. (2003) with both the Padova (Girardi et al. 2000) and Yale- 
Yonsei isochrones (Yi et ah 2001; Demarque et ah 2003). Previous similar studies allowed 
to confine the cluster age in the relatively large interval between 8 and 12 Gyrs (Carraro et 
ah 1994, Chaboyer et ah 1999, King et ah 2005, Carney et ah 2005, and reference therein). 
As discussed in Stetson et al (2003), and confirmed in this work, the cluster reddening is 
E(B-V) = 0.09±0.04. There is a large scatter in the literature on the absolute distance 
modulus (m — M) a estimates, which range in the interval between 12.6 and 13.6. These 
large uncertainties in both age and distance have been usually ascribed to uncertainties in 
the cluster metal abundance. The new specroscopic data presented in this paper allow to 
put on a more solid basis these fundamental parameters. 


6.1. Padova Isochrones 

Our [Fe/H] = 0.39 empirical measurement, translates into a metallicity Z = 0.046 (Car¬ 
raro et ah 1999) and implies a AY/AZ close to 2. We generated isochrones for this metal¬ 
licity, and for ages ranging from 7 to 11 Gyrs from Girardi et ah (2000). An appropriate and 
meaningful isochrone fit implies that all the loci of the CMD, e.g. the turn-off point (TO), 
the sub giant branch (SGB), the red giant branch (RGB), and the clump of He-burning stars 
must be simultaneously overlapped by the models. Our best fit (by eye) estimate is shown 
in Fig. 5 and 6, both in the V vs (B-V) and V vs (V-I) plane. In Fig. 5 we plot the V vs 
(B-V) CMD of NGC 6791, and superpose the whole set of isochrones, whereas, in Fig. 6, we 
only show the best fit isochrone in the V vs (B-V) and V vs (V-I) plane. 

The isochrone solutions in Fig. 5 have been obtained by shifting the theoretical lines by 
E(B — V) = 0.09 and (m — M)v = 13.35. Clearly ages older than 9 Gyr can be ruled out, 
since a fit to the TO with an older isochrones implies to decrease the distance modulus, but, 
in this way, the theoretical clump would be brighter than the observed one. On the other 
hand, also ages younger than 8 Gyr do not seem possible, since a fit of the TO region with a 
younger isochrone would result in a RGB redder than the observed one (impling a reddening 
value significantly larger that the observational limits), and also the clump magnitude would 
be fainter than the observed counterpart. 

Only the isochrones for ages of 8 and 9 Gyr provide a good fit. In details, the 9 Gyr 
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isochrone fits well the CMD with the adopted parameters, although the clump luminosity 
turns out to be slightly brighter than the observed one. On the other hand, the 8 Gyr 
isochrone must be shifted by E(B — V) — 0.09 and (m — M)y = 13.45 to provide a very 
good fit. This is shown in Fig. 6. We note that the lower MS is mismatched. This is a well 
known problem for metal rich clusters, as extensively discussed by Bedin et al. (2001), and 
it is likely dne to problems in the transformation of the models from the theoretical to the 
observational plane. 

On the overall, however, the fit is very good, and implies for NGC 6791 this set of fun¬ 
damental parameters: 8.0±1.0 Gyr, 13.07±0.05 (internal error), and 0.09±0.01 for the age, 
absolute distance modulus, and reddening, respectively. The associated errors are internal 
errors and have been estimated by eye. They simply reflect the degree of freedom we have 
to displace the isochrones, still achieving an acceptable fit. 


6.2. Yale-Yonsei Isochrones 

An independent determination of the age, distance, and a constraint on the reddening 
of NGC 6791 can be derived using the Y' 2 isochrones (Yi et al. 2001; Demarque et al. 
2003). The Padova and Y 2 isochrones were both constructed using the same OPAL opacities 
tables. Otherwise, the description of the microscopic and macroscopic physics, as well as 
the numerical procedures, differ in many details in the two sets of isochrones. The color 
transformations are also independently derived. 

In the Y 2 system, in which Z e = 0.0181, [Fe/H] = 0.39 corresponds to (Y, Z) = (0.31, 
0.04). The fit is based on the main sequence (MS) position just below the TO, the position 
of the TO point, the SGB, and the RGB color. A good fit, as shown in Fig. 7, is obtained 
for (m — M)y = 13.35, 0.1 mag smaller than the distance modulus we used for the Padova 
isochrone fit. The best fit is obtained by assuming a reddening E(B — V)= 0.13, somewhat 
larger than the reddening adopted in the previous section, but still within the estimated 
range of previous investigations. The age we derive from the Y 2 fit is between 8 and 9 
Gyr, in good agreement with the Padova age, even though one notes that the position of 
the lower main sequence differs in the two sets of isochrones. The unevolved main sequence 
of the Padova isochrones has a steeper downward slope than the observations, whereas the 
opposite holds for the Y 2 isochrones. 

Similarly, Chaboyer et al. (1999) conclude that the cluster age is 8.0 ± 0.5 Gyr, assuming 
[Fe/H] = 0.4, but using an older observational data set (Kaluzny & Rucinski 1995), and a 
version of the Yale stellar evolution code that slightly differs from the one used to construct 
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the T' 2 isochrones. In an analysis of their infrared photometry, Carney et ah (2005) derive 
an age between 9 Gyr (for [Fe/H] = 0.3) and 7.5 Gyr (for [Fe/H] = 0.5), also in good 
agreement with onr result. Both the Chaboyer et al. (1999) and the Carney et al. (2005) 
ages are consistent with the Padova and Y 2 fits described in this paper. We should note 
however that the Carney age estimates were obtained using the same set of Y 2 isochrones 
that we used in the present work, and therefore their age determination is not completely 
independent form our one. 

Stetson et al. (2003) derive a much older age (12 Gyr) with the help of unpublished 
VandenBerg isochrones. It appears that the large difference in age is due in part to the au¬ 
thors’ choice of a markedly smaller absolute distance modulus (12.79 mag). A superposition 
of the F' 2 isochrones for the range 8-12 Gyr, shifted by (m — M) 0 = 12.79 and E(B — V) 
= 0.09 (Stetson et al. adopted values), on the CMD of NGC 6791 is shown in Fig. 8, for 
comparison purpose. Although it is not possible to rule out completely the Stetson et al. 
(2003) fit, the disagreement with other well calibrated isochrones raises questions about the 
calibration of the VandenBerg isochrones. 

An additional, independent age and distance estimate is in King et al. (2005) who 
obtained an excellent fit of the upper main sequence, TO, and SGB both of Stetson et al. 
(2003) groundbased CMD of NGC 6791, and of their ACS/HST CMD in the F606W, F814W 
bands by using the Teramo isochrone set by Pietrinferni et al. (2004). From both fits, King 
et al. (2005) derived an age of 9 ± 1 Gyr, an absolute distance modulus (m — M) 0 = 13.0, 
and a reddening E(B-V)=0.12 for a metallicity [M/H]=+0.4, and Y=0.288. Also in the fit of 
the rriFsuw vs. nTF 606 W-^F 8 i 4 W ACS/HST CMD, the Teramo isochrones tends to be redder 
and redder going to fainter magnitudes, starting from ~ 2 magnitudes below the TO, as 
already noticed for the Padova isochrones. 

Finally, we must mention that the precise age of NGC 6791 also depends on the adopted 
value of the ratio (AY/AZ) for galactic helium enrichment. This quantity is poorly known; 
it may be a function of Z, and may differ from system to system. Demarque et al. (1992) 
have found that varying Y from 0.32 to 0.36 could reduce the age of the cluster by as much 
as 15%. The age estimate of NGC 6791 might have to be increased if the enrichment ratio 
AY/AZ is much less than 2 (isochrones for AY/AZ near 2 were assumed in the Padova, 
Y 2 , Teramo, and VandenBerg fits). 

With our present knowledge of the AY/AZ parameter, we conclude that the age of 
NGC 6791 must be in of range 7.5-8.5 Gyr, with a higher preference towards the higher 
limits. We point out that, even in the unlikely event that the age of NGC 6791 is as low 
as 7.5 Gyr, its high metallicity and age present a major challenge to the accepted view of 
Galactic chemical enrichment. 
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In conclusion, three sets of independent isochrones consistently inply that the age of NGC 
6791 is around 8 Gyr, adopting the metallicity and the reddening coming from observa¬ 
tions. The difference in reddening might simply ascribed to the different Helium abundance 
adopted, and to some photometric zero point error. 


7. DISCUSSION and CONCLUSIONS 

The hardest point with NGC 6791 is how inside this cluster such high metallicity stars 
could have been produced. In fact, this cluster does not have any counterpart in the Milky 
Way. We note here that Kinman (1965) originally identified NGC 6791 as a globular clus¬ 
ter. Even if this interpretation were to be adopted, the high metallicity of NGC 6791, much 
higher than that of any Galactic globular cluster or nearby dwarf galaxy, remains mysterious. 


With Berkeley 17 and Collinder 261, NGC 6791 is one of the oldest open clusters of the 
Galaxy (Carraro et al. 1999), but its metal abundance is incomparably higher. Moroever, 
NGC 6791 is one of the most massive open clusters (4000 M e at least). It lies at 1 kpc 
above the Galactic plane, inside the solar ring. This combination of mass and position is 
hard to explain, since the interaction with the dense Galactic environment should strongly 
depopulate a typical open cluster. 


NGC 6791 is routinely considered in the studies of the chemical evolution of the Galac¬ 
tic disk, and occupies a unique position in the Galactic disk radial abundance gradient (see 
Fig. 9). By including NGC 6791, the slope of the gradient changes from -0.05 (solid line) to 
-0.07 (dashed line). Besides, if one considers the slope defined only by clusters older than 4 
Gyr (Friel et al. 2002, Fig. 3, upper panel), the slope doubles, from -0.06 to -0.11. 

In the same figure, the horizontal solid line indicates the epicyclical amplitude of NGC 6791 
orbit (see below, and Carraro & Chiosi 1994). One can readily see how NGC 6791 is quite 
an exotic object. If, by chance, at the present time the cluster would be at different orbit 
phase, which would put it, e.g., beyond 12 kpc, there would be a drastic change and even an 
inversion of the slope of the Galactic disk abundance gradient. Finally, the position of this 
cluster in the Galactic disk Age-Metallicity relationship (Carraro et al. 1998) is puzzling as 
well, since the cluster significantly deviates from the mean trend. 


In Fig 10 we present NGC 6791 Galactic orbit. This was obtained by integrating back 
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in time ( 1 Gyr) the cluster from its present position and kinematics using the Galaxy 
N-body/gasdynamical model by Fux (1997, 1999). The adopted radial velocity and proper 
motions come from Geisler (1988) and Cudworth (private comunication), whereas the Galac- 
tocentric rectangular initial conditions (positions and velocities) was derived as in Carraro 
& Chiosi (1994). 

Intererstingly enough, this plot shows that the cluster moves from the outer disk regions of 
the Milky Way, more than 20 kpc far away from the Galactic center, and enters the Solar 
Ring going as close as 6 kpc from the Galactic center. 


The eccentricity (e=0.59) of this orbit is quite high for a Population I star cluster 
(Carraro & Chiosi 1994), and it is much more similar to a globular cluster/dwarf galaxy 
orbit. 

A plausible scenario is that NGC 6791 is what remains (the nucleus) of a much larger 
system, which underwent strong tidal disruption. This would explain the cluster orbit, and 
provide a reasonable explanation for the high metallicity of its stars, which could have been 
produced only inside a deep potential well. 

However, within the observational errors, we did not find any significant abundance 
spread. This would mean that the bulk of the stars in the cluster was produced in a sin¬ 
gle burst of star formation. This fact makes more difficult the capture interpretation since 
Local Group Galaxies normally exhibit spreads in metal content and possess lower metal 
abundance (Mateo 1998). We stress however the fact that our results are based on only ten 
stars, and that only larger spectroscopic surveys can better address this particular problem. 


An alternative more conservative scenario is that the cluster was born in the inner side 
of the Galaxy, close to the bulge, where the metal enrichment has been fast. Grenon (1999) 
studied the kinematics of a group of old (10 Gyrs) metal rich [M/H] > 0.30 stars and 
suggested that they formed close to the bulge and then migrated at large Galactocentric 
distance due to the perturbation of the Galactic bar. 

The orbit we calculated actually includs the effect of the bar, and NGC 6791 indeed 
moves well outside the solar circle. NGC 6791 is very concentrated for an open cluster, and 
spent most of its time at moderate Galactic latitude. This might help to explain its survival. 


The observations described in this paper were carried out remotely from Yale University 
by Giovanni Carraro. We deeply acknowledge Diane Harmer, George Will and Chris Hunter 
for support and help. The work of GC is supported by Fundacion Andes. GP and SV 
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Fig. 5.— Five Padova isochrones for the ages of 7, 8, 9, 10, and 11 Gyr are fitted to the 
observed CMD in the V vs (B-V) plane. The isochrones are shifted by E(B — V) — 0.09 
and (m — M) v = 13.35. 
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(B-V) (V—I) 


Fig. 6.— Best fit isochrone solution of the CMD of NGC 6791 with the Padova models. The 
isochrone and setting parameters are indicated in the plot. 
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(B-V) 


Fig. 7.— Yale-Yonsei isochrone solution for ages of 8 and 9 Gyr. 
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(B-V) 


Fig. 8.— Yale-Yonsei isochrone solution for ages of 8, 10 and 12 Gyr (from the top to the 
bottom). Note how also the two older age isochrones are clealy ruled out. 
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^(GC) 


Fig. 9.— The Galactic disk chemical abundance radial gradient. The data are from Fricl 
et al. (2002), with the exception of Berkeley 22 and Berkeley 66 taken from Villanova et al. 
(2005) and NGC 6791 (filled square, coming from the present study). The solid line is the 
linear fit without NGC 6791, whereas the dashed line is a linear fit to all the data points. 
The horizontal dotted line shows the epicyclical amplitude of NGC 6791 orbit. 
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Fig. 10.— Galactic orbit of NGC 6791 in the X-Y and meridional plane. The position of 
the Sun is indicated. 
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Table 1. OBSERVED STARS 


ID 

K65 

RA 

DEC 

V 

(B-I) 

Vrad (km s b 

S/N 

Spectral Type 

10898 


19:21:01.13 

+37:42:13.80 

14.459 

3.059 

-47.37+0.088 

40 

K4 III 

11814 

3003 

19:21:04.27 

+37:47:18.90 

13.849 

3.233 

-46.52+0.082 

60 

K4/5 III 

1442 


19:21:16.33 

+37:52:15.80 

14.056 

2.580 

-11.84+0.056 

100 


2044 


19:20:30.92 

+37:48:45.30 

14.144 

2.917 

-5.27+0.072 

40 


2423 


19:20:33.26 

+37:50:12.70 

14.082 

2.794 

-220.41+0.066 

70 


2793 

3036 

19:20:34.86 

+37:46:30.10 

14.538 

2.751 

-48.00+0.108 

30 

Ml/2 III 

3369 

3030 

19:20:37.89 

+37:44:49.30 

14.529 

2.673 

-48.37+0.094 

40 

K2/3 III 

4162 


19:20:40.85 

+37:46:21.80 

14.551 

2.687 

-52.11+0.146 

10 

K2/3 III 

4715 


19:20:42.73 

+37:51:07.70 

14.515 

2.698 

-47.91+0.088 

50 

K2/3 III 

5583 


19:20:45.58 

+37:39:51.20 

14.602 

2.742 

-43.62+0.196 

10 

K3/4 III 

6940 

3013 

19:20:49.67 

+37:44:08.00 

14.588 

2.659 

-46.00+0.155 

10 

K2/3 III 

7922 


19:20:52.47 

+37:50:15.80 

14.482 

2.671 

-48.28+0.080 

30 

K2/3 III 

7972 

3010 

19:20:52.60 

+37:44:28.50 

14.136 

3.356 

-44.35+0.094 

40 

K7 III 

8082 

SE-49 

19:20:52.89 

+37:45:33.40 

14.546 

2.639 

-46.18+0.102 

40 

K2/3 III 

8266 

2001 

19:20:53.39 

+37:48:28.40 

13.741 

3.395 

-47.73+0.080 

40 

K9 III 

852 


19:20:22.40 

+37:51:42.40 

14.738 

2.748 

-67.91+0.267 

10 


8563 


19:20:54.19 

+37:46:28.80 

14.554 

3.071 

-42.44+0.084 

40 

K4 III 

8904 

2008 

19:20:55.11 

+37:47:16.50 

13.862 

3.603 

-46.91+0.098 

30 

MO III 

8988 

3018 

19:20:55.31 

+37:43:15.60 

14.557 

3.005 

-47.46+0.088 

30 

K4 III 

95 


19:20:11.19 

+37:49:48.70 

13.589 

2.986 

-72.93+0.076 

50 
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Table 2. ADOPTED ATMOSPHERIC PARAMETERS 








Table 3. LINELIST 


6408.016 Fel 
6439.075 Cal 
6469.123 Fcl 
6491.561 Ti2 
6499.650 Cal 

6554.230 Til 
6575.020 Fel 

6609.120 Fel 
6634.100 Fel 
6677.955 Fel 
6705.131 Fel 
6721.848 Sil 

6743.120 Til 


6411.650 Fel 
6449.808 Cal 
6469.210 Fel 

6494.980 Fel 
6501.691 Til 
6556.070 Til 
6581.221 Fel 
6625.041 Fel 
6643.640 Nil 
6677.990 Fel 
6710.310 Fel 
6725.390 Fel 
6743.185 Til 


6419.980 Fel 
6455.598 Cal 
6471.662 Cal 
6496.897 Ball 
6518.380 Fel 

6569.230 Fel 
6593.880 Fel 
6627.560 Fel 

6646.980 Fel 
6698.673 All 
6713.760 Fel 
6733.160 Fel 
6750.150 Fel 


6421.350 Fel 


6527.202 Sil 
6572.779 Cal 
6606.970 Till 
6633.440 Fel 

6663.231 Fel 
6703.570 Fel 
6715.410 Fel 

6737.980 Fel 


6436.430 Fel 
6464.661 Fel 
6481.880 Fel 
6498.950 Fel 
Nil 

6574.240 Fel 
6608.030 Fel 
6633.760 Fel 
6663.450 Fel 
6705.101 Fel 
6717.681 Cal 
6741.628 Sil 





Table 4. 


MEAN STELLAR ABUNDANCES 


ID 

[M/H] 

[Fe/H] 

[Cal/H] 

[Ti/H] 

[Ba/H] 

[Si/H] 

[Ni/H] 

[Al/H] 

10898 

0.37±0.06 

0.38d=0.08 

0.33=1=0.12 

0.38=1=0.02 


0.42 

0.32 

0.20 

11814 

0.39+0.09 

0.34d=0.08 

0.32=1=0.08 

0.34=1=0.05 


0.39=1=0.03 

0.29=1=0.14 

0.23 

3369 

0.41±0.09 

0.37±0.09 

0.35=1=0.05 

0.38=1=0.07 

0.19 

0.42 

0.36 

0.18 

4715 

0.39=1=0.07 

0.36±0.05 

0.33=1=0.07 

0.35=1=0.02 

0.18 

0.41=1=0.02 

0.42 

0.23 

7922 

0.39=1=0.02 

0.39±0.07 

0.33=1=0.14 

0.36=1=0.02 


0.37=1=0.09 

0.40 

0.22 

7972 

0.40=1=0.09 

0.39±0.05 

0.35=1=0.06 

0.34=1=0.05 

0.20 

0.39 

0.41 

0.26 

8082 

0.42=1=0.04 

0.38+0.05 

0.38=1=0.07 

0.38=1=0.07 

0.14 

0.36 

0.48 

0.23 

8266 

0.37=1=0.05 

0.37+0.05 

0.38=1=0.07 

0.37 

0.24 

0.39 

0.31 

0.19 

8563 

0.38=1=0.08 

0.38=1=0.06 

0.32=1=0.03 

0.36=1=0.04 

0.32 

0.42 

0.37 

0.25 

8988 

0.36=1=0.06 

0.40=1=0.07 

0.35=1=0.10 

0.31=1=0.07 

0.22 

0.39 

0.34 



Table 5. ABUNDANCE RATIOS 


ID 

[Fe/H] 

[Ca/Fe] 

[Ti/Fe] 

[Ba/Fe] 

[Si/Fe] 

[Ni/Fe] 

[Al/Fe] 

10898 

0.38=1=0.08 


0.00 


+0.04 

-0.06 

-0.18 

11814 

0.34=1=0.08 

-0.02 

0.00 


+0.05 

-0.05 

-0.11 

3369 

0.37=1=0.09 

-0.02 

+0.01 

-0.14 

+0.05 

-0.01 

-0.19 

4715 

0.36=1=0.05 

-0.03 

-0.01 

-0.15 

+0.05 

+0.06 

-0.13 

7922 

0.39=1=0.07 

-0.06 

-0.03 


-0.02 

+0.01 

-0.17 

7972 

0.39=1=0.05 

-0.04 

-0.05 

-0.15 

0.00 

+0.02 

-0.13 

8082 

0.38=1=0.05 

0.00 

0.00 

-0.22 

-0.02 

-0.02 

-0.15 

8266 

0.37=1=0.05 

+0.01 

0.00 

-0.12 

+0.02 

+0.02 

-0.18 

8563 

0.38=1=0.06 

-0.06 

-0.06 

0.00 

+0.04 

-0.01 

-0.13 

8988 

0.40=1=0.07 


-0.06 

-0.13 

-0.01 

-0.01 














